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1. Summary 

Progression through the eukaryotic cell cycle is characterized by specific tran- 
sitions, where cells move irreversibly from stage i — 1 of the cycle into stage i. 
These irreversible cell cycle transitions are regulated by underlying bistable 
switches, which share some common features. An inhibitory protein stalls pro- 
gression, and an activatory protein promotes progression. The inhibitor and 
activator are locked in a double-negative feedback loop, creating a one-way 
toggle switch that guarantees an irreversible commitment to move forward 
through the cell cycle, and it opposes regression from stage i to stage i - 1. In 
many cases, the activator is an enzyme that modifies the inhibitor in multiple 
steps, whereas the hypo-modified inhibitor binds strongly to the activator 
and resists its enzymatic activity. These interactions are the basis of a reaction 
motif that provides a simple and generic account of many characteristic prop- 
erties of cell cycle transitions. To demonstrate this assertion, we apply the 
motif in detail to the Gl/S transition in budding yeast and to the mitotic check- 
point in mammalian cells. Variations of the motif might support irreversible 
cellular decision-making in other contexts. 

2. Introduction 

During the cell division cycle, a cell's DNA molecules are accurately replicated 
and then each pair of identical DNA molecules are carefully partitioned to 
the pair of newly forming daughter cells, so that each daughter cell gets one 
and only one copy of each DNA molecule. The perpetuation of life demands 
that these processes be carried out with nearly flawless precision, although 
some errors are unavoidable, of course. In addition to the complex biochemi- 
cal and mechanical mechanisms that carry out DNA replication and cell 
division, cells use subtle biochemical regulatory networks to ensure the 
proper completion and sequencing of cell cycle events. 

In eukaryotic cells, the cell cycle regulatory network is based on transcrip- 
tional and post-translational controls of a class of cyclin-dependent protein 
kinases (CDKs). As cells proceed periodically through the classic phases of 
the cell cycle (Gl, S, G2, M, Gl, . . .), these CDKs and their associated regulators 
are periodically activated and inactivated, thereby driving progression from one 
stage of the cycle to the next. Proper progression through the cell cycle requires 
that stage i must not begin until stage i — 1 is successfully completed and that 
the transition from stage i — 1 to stage i is irreversible. 

In general, there are five major transitions at boundaries between sequential 
stages of the cell cycle in eukaryotes: Gla/Glb (the 'restriction point 7 in mam- 
malian cells, 'Start' in yeast), Gl/S (initiation of DNA synthesis), G2/M (entry 
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Figure 1. Generic picture of a bistable cell cycle switch, (a) Signal -response curve. The activator of stage / of the cell cycle ('response') is a bistable function of a 
'signal' generated by stage / - 1, for 0 < signal < threshold. When the signal exceeds the threshold, the activator turns on and the cell proceeds to stage /'. 
Subsequently, even if the signal drops to zero, the activator stays on and the cell remains in stage /'. (b) An activator- inhibitor pair involved in a double-negative 
feedback loop can generate the sort of bistable response postulated in (a). (c,d) Phase planes for an activator- inhibitor system for (c) signal < threshold and for 
(d) signal > threshold. 



Table 1. Activators and inhibitors of certain cell cycle transitions. S, stoichiometric inhibitor; E, enzymatic inhibitor. 



transition 


activator 


inhibitor 


restriction point 


E2F transcription factor 


retinoblastoma protein (S) 


Start (yeast) 


SBF transcription factor 


Whi5 (S) 


G1/S (yeast) 


Clb5-dependent kinase 


Sid (S) 


G1/S (yeast) 


Clb5-dependent kinase 


Cdh1 (E) 


G2/M 


cyclin B-dependent kinase 


Wee1 (E) 


metaphase/anaphase 


anaphase-promoting complex 


mitotic checkpoint complex (S) 



into mitosis), metaphase/anaphase (chromosome alignment 
on the mitotic spindle) and telophase /Gl (exit from mitosis 
and return to Gl phase). Progression through each of these 
transitions is controlled by a balance between inhibitors 
and activators (table 1) [1]. Before the transition, the inhibitor 
overwhelms the activator. As soon as the prior stage is 
judged to be complete, a signal helps the activator to elimin- 
ate the inhibitor and the transition occurs. The transitions are 
irreversible because each biochemical mechanism functions 
as a 'bistable' switch (figure la). The connection between irre- 
versible cell cycle transitions and bistable biochemical 
switches has been emphasized repeatedly in the theoretical 
literature [2-5] and has been confirmed experimentally in 
all cases [6-11]. 

In each case, the bistable switch could potentially provide 
'checkpoint' control over the cell cycle transition. By monitor- 
ing defective states of chromosomes (damaged, not fully 
replicated, not fully aligned on the mitotic spindle), intra- 
cellular surveillance mechanisms could keep the signal 
(figure la) at a subthreshold level, thereby blocking the 



transition. When the defect is repaired, the signal is released 
to drive the transition. In any given cell type, only a few of 
the five cell cycle transitions function as 'checkpoints' in 
this sense. But each checkpoint mechanism (we propose) 
relies on the bistability of the underlying transition. 

Bistability at cell cycle transitions is the consequence of 
a double-negative feedback loop by which the activator 
inactivates its own inhibitor (figure lb). To illustrate the con- 
nection between double-negative feedback loops and 
bistability /irreversibility, we provide two generic 'phase - 
plane diagrams' in figure lc,d (see the electronic supplemen- 
tary material, §S4 for a description of how to interpret phase 
plane diagrams). When signal = 0 (figure lc), there are two 
stable steady states of the control system: 'checkpoint engaged' 
(lots of inhibitor, very little functional activator), and 'check- 
point disengaged' (very little inhibitor, lots of activator). At 
the beginning of the process, the signal is off and the check- 
point is engaged. As the signal increases from zero towards 
the threshold, the level of functional activator increases only 
a little (figure la), and the checkpoint remains engaged. As 



soon as the signal exceeds the threshold, the lower steady state 
(checkpoint engaged) disappears, and the control system is 
forced to transition to the upper steady state (checkpoint 
disengaged), as in figure Id. After the transition, even if the 
signal drops to zero (i.e. return to figure lc), the checkpoint 
mechanism remains locked in the stable, disengaged steady 
state. In this sense, the transition is irreversible. 

Each of the five cell cycle transitions mentioned previously 
is controlled by its own unique biochemical mechanism. In 
some cases, the inhibitor is an enzyme. For example, at the 
G2/M transition, the protein kinase Weel phosphorylates 
and inhibits the mitotic M-CDK. During the transition, Weel 
gets multiply phosphorylated and inactivated by the active 
(dephosphorylated) form of the M-CDK [12]. 

In other cases, a stoichiometric inhibitor binds to the activa- 
tor. For example, before the restriction point in mammalian 
cells, the retinoblastoma protein (RB) binds to and inhibits 
members of the E2F family of transcription factors [13]. As 
cells pass the restriction point, some E2F molecules escape 
inhibition by RB and begin to activate expression of the 
CYCE gene. The rise of cyclin E-dependent kinase hyper- 
phosphorylates/ inactivates RB, thereby instantiating the 
regulatory motif in figure lb. The same is true at Start in 
yeast. Before Start, Whi5 binds to and inhibits the transcription 
factor, SBF. At the transition point, SBF turns on production of 
the Gl cyclins, Clnl and Cln2, and the Cln-dependent kinases 
phosphorylate and neutralize Whi5 [14,15]. 

The Gl/S transition is also controlled by a stoichiometric 
inhibitor: p27 Kipl in mammalian cells, p40 Slcl in budding 
yeast and p25 Ruml in fission yeast [16]; generically known as 
cyclin-dependent kinase inhibitors (CKIs). These CKIs bind 
to and inhibit S phase-specific CDKs (cyclin A:Cdk2 in mam- 
malian cells, Clb5:Cdkl in budding yeast and Cig2:Cdkl in 
fission yeast). At the Gl/S transition, CKI is phosphorylated 
and degraded, releasing the S-CDK to initiate DNA synthesis. 
According to our paradigm (figure 1 and table 1), bistability at 
the Gl/S transition requires that CKI degradation be promoted 
by the S-CDK complex. CKI degradation is indeed pro- 
moted by CDK-mediated hyper-phosphorylation [17-19]. 
The abrupt degradation of Sicl at the Gl/S transition in 
budding yeast is usually attributed to multi-site phosphoryl- 
ation of Sicl by Cln2-dependent kinase, which is resistant to 
Sicl inhibition [20]. However, given that Cln2 is degraded in 
S phase, phosphorylation of Sicl by Cln2-dependent kinase 
alone cannot account for irreversibility of the transition. 
Irreversibility of the Gl/S transition requires that Sicl be 
hyper-phosphorylated by Clb-dependent kinases as well [21], 
a supposition recently confirmed by Koivomagi et al. [22]. 

Finally, the meta /anaphase transition is induced by the 
anaphase-promoting complex /cyclosome (APC/C). Before 
the transition, the APC is kept inactive by binding with its stoi- 
chiometric inhibitor, the mitotic checkpoint complex (MCC). 
(The MCC is a large protein complex, including Mad2 and 
Cdc20 among other components.) Bistability of the checkpoint 
mechanism requires that the MCC is inactivated by the APC. 
The biochemical basis of MCC inactivation during anaphase 
was enigmatic for a long time, until Reddy et al. [23] showed 
that disassembly of the MCC: APC complex is promoted 
by poly-ubiquitination of the Cdc20 subunit of the MCC, 
catalysed by the APC itself. Although these findings were 
questioned later [24], more recent publications have confirmed 
the conclusions of Kirschner's laboratory [25-27]. Briefly, APC 
promotes poly-ubiquitination of the Cdc20 molecule in the 



MCC complex, which leads to Cdc20 degradation, MCC disas- 
sembly and release of free APC that associates with another 
MCC or Cdc20. 

Hence, in both cases, CKI and MCC are substrates as well 
as inhibitors of the very enzyme complexes they are opposing 
(figure lb), thereby creating the potential for bistability of the 
molecular mechanisms and irreversibility of the associated 
cell cycle transitions. Furthermore, the enzymes carry out 
multiple modifications of their substrate /inhibitor (multi- 
site phosphorylation of CKIs and poly-ubiquitination of the 
Cdc20 subunit of the MCC). 

Based on the experimental evidence summarized here, 
we present in this paper a general and robust mechanism of 
cell cycle transitions based on bistable biochemical switches. 
The mechanism relies on multiple modifications of an inhibi- 
tor by the enzyme it inhibits. The inhibitor strongly inhibits 
its first modification by the enzyme, but later it becomes a 
good substrate for further modifications. We characterize the 
properties of this mechanism by mathematical modelling, 
we describe a number of realizations of this mechanism in 
cell cycle controls, and we compare the model's properties 
with the observed genetic and physiological properties of 
cell cycle transitions. 



3. Results 

3.1. Substrate -inhibitor multiply modified: a general 
motif for bistability at cell cycle transitions 

We start with a simple model (figure 2a) of an inhibitor I that 
is multiply modified by an enzyme A, which is an activator of 
a cell cycle transition. IM represents 'early' modifications 
of the inhibitor, which are sluggishly accomplished and 
which tie up the activator in tightly bound enzyme -substrate 
complexes, A:I. IMM represents 'later' modifications, which 
occur quickly and (perhaps) processively. It is important 
that the enzyme disengage from its substrate at some point 
between the early and later modifications. This distribu- 
tive step [28] is a consequence of our assumption that A:I 
is tightly bound and slowly processed, whereas A:IM is 
loosely bound and rapidly processed. (The model is bistable 
as described here even if as much as approx. 20% of the 
dual modification process (I — > IM — > IMM) occurs proces- 
sively rather than distributively, as shown in the electronic 
supplementary material, §S1.) 

We assume that the enzymatic modifications follow the 
Michaelis-Menten mechanism. Because enzyme A and sub- 
strate I are present in comparable concentrations, we must 
implement the Michaelis-Menten mechanism in terms of 
elementary monomolecular and bimolecular steps, obtaining 
a set of nonlinear ordinary differential equations (ODEs) 
in terms of I = [I], I M = [IM], I MM = [IMM], C = [A:I], and 
C M =[A:IM]: 



dl 

dt ~ 
dlu 
dt 

dluN 
~df 



-k asl A-I + k dil C + k dml h 



-kaslA • Im + fcd&QM + kcatlC ~ ^dml^M + ^dm2-5v 



: ^cat2CM — fcdm2^MM, 



(3.1) 

(3.2) 
(3.3) 
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(c) (d) 
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Figure 2. The SIMM motif, (a) Wiring diagram. Enzyme A binds to substrate I to form an enzyme -substrate complex, A:l, which then reacts to form a modified 
substrate IM plus free enzyme. Then, the same sequence of steps occurs a second time to form the doubly modified form of I. (b) Signal -response curve. For a 
fixed value of total enzyme concentration, Aj = 1, the steady-state concentration of free enzyme, A SSI is plotted as a function of total substrate concentration, / T . 
Path b is explained in (rf). (c) Signal -response curve. For a fixed value of total substrate concentration, / T = 4, the steady-state concentration of free enzyme, A SSI 
is plotted as a function of total enzyme concentration, 4 T . Path c is explained in (d). (d) Control plane. As a function of 4 T and / T , the reaction network is bistable 
only within the V-shaped region. Below the V, the checkpoint is engaged (A inhibited), and above the V, the checkpoint is disengaged (A active). The black circle 
represents the neutral state of the checkpoint before the cell cycle transition: the reaction system is in the bistable region, in the checkpoint-engaged stable 
steady state. Paths b and c are possible paths for disengaging the checkpoint and triggering an irreversible cell cycle transition, as shown in (b,c). Indeed, any 
path that carries the control system across the upper boundary of the V and back to the neutral state will trigger an irreversible cell cycle transition. Eventually, 
the checkpoint must be reset to the 'engaged' state, which requires the system to follow a path r that crosses the lower boundary of the V and then returns to 
the neutral state. 



— = k asl A ■ I - k dil C - k catl C, (3.4) 

— jj^ = ^as2^ ' — fcd&QM — ^cat2Qy[ (3.5) 

dA 

and — = -k as iA • I + k dil C + k catl C - k as2 A ■ I M 

+ kdilCyi + k cat2 C M - (3.6) 

Subscripts on rate constants have the following meanings: 
as, association; di, dissociation; cat, catalysis; dm, de- 
modification; 1, first modification; 2, second modification. 
Note that equations (3.1)- (3.6) satisfy two conservation con- 
ditions: A T = A + C + C M = const, and J T = I + C + Im + 
Cm + ^mm = const. We study the behaviour of this reaction 
motif for the rate constant values given in table 2. These 
values are consistent with our assumptions that I binds 
very strongly to A to form a stable A:I complex, which 
reacts slowly to form the early-modified forms IM, and that 
the later-modified forms are good substrates — not good 
inhibitors — of A. 

As simple as this mechanism appears, it is a double-negative 
feedback loop (figure lb). As long as A T is sufficiently less than J T , 
most of the activator is locked up in the inactive A:I complex. 
However, when A T > I T , there will be some free enzyme 
(A = A T — J T ) that can catalyse the later modification step, 



IM — > IMM. Because the later-modified forms of I are not good 
inhibitors of A, the proportion of functional activator in the reac- 
tion network will increase steadily. Indeed, as shown by the 
signal -response curves (figure 2b/:), the network motif exhibits 
bistability for the rate constant values proposed in table 2a. 

We assume that, when a cell is in stage i — 1 of the cell 
cycle, J T is sufficiently greater than A T so that the control 
system is locked in the state of low activity of A (check- 
point-engaged). To release the checkpoint, the mechanism 
must generate a signal that either decreases J T or increases Aj. 
When the A T /I T ratio gets large enough, the checkpoint dis- 
engages, and the cell transitions to stage i of the cell cycle. 
(In figure 2b, J T decreases at fixed A T = 1, and the transition 
occurs at J T ~ 2.5; in figure 2c, A T increases at fixed J T = 4, 
and the transition occurs at A T w 1.6.) 

The bistable behaviour of this reaction mechanism depends 
critically on the total concentrations of the inhibitor/ substrate I 
and activator/ enzyme A, as shown in figure 2d. Bistability 
is restricted to a V-shaped region in the A T , J T plane (called a 
'control plane' because the parameters A T and J T control the 
behaviour of the reaction mechanism). Below the V, the 
system is mono-stable with little active A (checkpoint 
engaged), and above the V, the system is mono-stable with 
lots of active A (checkpoint disengaged). 

Bistability of the reaction motif in figure 2a depends on 
three properties of the checkpoint protein I: it is a substrate 



Table 2. Values of rate constants and other parameters in the models. U, concentration unit (approx. 10 nM); T, time unit (approx. 1 min). 



rate constants 


units 


interpretation 


(a) SIMM motif (figure 2) 


* as i = 100, k as2 = 50 


u" 1 r 1 


1 binds more strongly to A than does IM 


* dn = 0.5, * di2 = 1 


r 1 


dissociation rate constant is same for A:l and A:IM 


*cat1 = 0.5, k cat2 = 50 


r 1 


IM is a much better substrate for modification by A than is 1 


^dml = 1/ ^dm2 = 1 


r 1 


de-modification rate constant is same for IM and IMM 


Aj, lj variable 


u 


total concentrations of A and 1 can be varied 


(b) SIMM* motif (figure 3). Same as set (a), plus 


* m = 50 


U-1 T -i 


rate constant for A + IM -> A + IMM 


4i = 0.01 


u 


Michaelis constant for first modification reactions 


(c) Sic-Clb mechanism (figure 4a) 


k s = 0.2 


ur 1 


rate of synthesis of Sid 


k d = 0.05 


r 1 


background degradation rate for all Sid forms 


k ash = 100 


U" 1 T" 1 


binding of Sid to Clb:Cdk1 


*du» = 0.01 


r 1 


dissociation rate constant for ClbiSic trimers 


* P jb = 20 


it 1 r 1 


rate constant for Clb-dependent phosphorylation of SicP to SicPP 
(followed by rapid degradation of SicPP) 


k 6p = 0.5 


r 1 


dephosphorylation rate constant 


4t_b = 0.5 


r 1 


rate constant for ClbrSic -> Clb + SicP 


[ClbT] variable 


u 


total concentration of Clb:Cdk1 dimers 


(d) Cln- Sic-Clb mechanism (figure 4c). Same as set (c), plus 


^as_n — 2 


ir 1 r 1 


binding of Sid to Cln:Cdk1 


^di_n = 2 


r 1 


dissociation rate constant for ClnrSic trimers 


* P _n = 0.1 


u- 1 r 1 


rate constant for phosphorylation of SicP by Cln-kinase 


^cat_n = 2 


r 1 


rate constant for CIn:Sic — > Cln + SicP 


[ClnT] 


u 


total concentrations of ClmCdkl dimers 


(e) MCC-APC mechanism (figure 5a) 


K m = 0.01 


u 


Michaelis constant APCMCC complex 


*u = 10 


U-i T -i 


rate constant for ubiquitination of MCCU by APC 


*du = 5 


r 1 


rate constant for ubiquitination of MCCU by APC 


/r cat = 0.5 


r 1 


rate constant for APCMCC -> MCCU + APC 


*a«T = 1 


r 1 


rate constant for activation of Mad2 by tensionless centromeres 


X tens variable 




fraction of centromeres under tension on mitotic spindle 


*ssec = 0.1 


ur 1 


rate constant for synthesis of securin 


4ec = 0.05 


r 1 


rate constant for degradation of securin 


^dsec,apc 0.5 


ir 1 r 1 


rate constant for APC-mediated degradation of securin 


(f) CycB- MCC-APC mechanism (figure 5d). Same as set (e) (except k a Nj), plus 


*scyc = 0.01 


ur 1 


rate constant for synthesis of cyclin B 


*dcyc = 0.01 


r 1 


rate constant for degradation of cyclin B 


^dcyc_apc 1 




rate constant for APC-mediated degradation of cyclin B 


*an = 1 


U-i T -i 


rate constant for activation of tensionless centromeres by CycB 


k in [CAPP] = 2 


r 1 


rate constant for inactivation of tensionless centromeres by CAPP 


*as = 100 


u- 1 r 1 


rate constant for association of APCMCC complex 


/fdi = 0.5 


r 1 


rate constant for dissociation of APCMCC complex 


kJh = 5 


r 1 


rate constant for activation of Mad2 by tensionless centromeres 



for enzymatic modification by A and an inhibitor of enzyme (SIMM) motif. We propose that SIMM is a general motif 

A as well, and I is multiply modified by A. Hence, we refer for cell cycle checkpoints with the characteristics described 

to figure 2a as a substrate -inhibitor-multiply modified in figure 1. 




Figure 3. The SIMM* motif, (a) Wiring diagram. For the second modification of I, we replace the Michaelis-Menten mechanism by a simple mass-action rate law. 
(b) Control plane, Aj versus / T . The V-shaped bistability zone for the SIMM* motif (solid line) is almost identical to the zone for the full SIMM motif (dashed line, 
from figure 2d), (c) Phase plane, / MM versus /, in the bistable zone (Aj = 2 and / T = 6). There are two stable steady states — checkpoint engaged (1) and check- 
point disengaged (3) — separated by an unstable steady state (2). (d) Phase plane in the monostable zone (Aj = 3 and / T = 6). There is now a single steady state: 
the stable checkpoint-disengaged state (3). In (c,d), the arrows indicate the direction of change of the state variables, / MM and I, as predicted by the dynamical 
system, equations (3.8) and (3.9). 



In figure la we stressed that a cell cycle checkpoint must be 
non-resettable (i.e. the signal that represents completion of 
stage i — 1 must trigger the checkpoint transition, and when 
the signal drops to zero after the transition, the checkpoint 
mechanism must remain in the disengaged state). It would 
appear from figure 7b jc that our SIMM motif does not have 
this property, but that depends on the starting values of J T 
and A T when the signal = 0. We suggest that, when the 
signal = 0, the 'neutral' values of J T and A T are in the bistable 
region (the black dot in figure 2d), and the control system 
is in the checkpoint-engaged steady state. As the signal 
increases, representing progression through stage i — 1, either 
I T decreases (path b) or A T increases (path c) until the check- 
point disengages (see the corresponding paths in figure 2b, c). 
If the signal drops to zero after the transition, the dynamical 
system returns to its starting point (•) on figure 2d, but the 
checkpoint is now locked in the disengaged steady state. 

Sometime later in the present cell cycle or early in the next 
cycle, the checkpoint must be reset to the engaged state (low 
activity of A). This requires a separate resetting mechanism, 
which takes the dynamical system along path r in figure 2d, 
from the bistable starting point (•) in the disengaged state 
(high activity of A) across the lower switching boundary and 
back to the starting point in the checkpoint-engaged state. 

The basic mechanism we are proposing for generating 
bistability by multiple modification of a checkpoint inhibi- 
tor was motivated by an earlier proposal by Kholodenko 
and co-workers [29,30] for bistability based on dual phos- 
phorylation of substrates by enzymes of the MAP kinase 
pathway. The reaction mechanisms (figure la; fig. 1 in [30]) 
are similar, but the details of the analysis are different. Mar- 
kevich et al. [30] assume that all enzymatic reactions are 



described by Michaelis-Menten kinetics, taking into account 
competition of the various substrate forms (S, SP and SPP) for 
the free kinase and phosphatase enzymes. Their use of 
Michaelis-Menten rate laws assumes that the three forms 
of the substrate are in great excess over the total amounts 
of kinase and phosphatase. Because the enzymes and sub- 
strates (A and I) operating at cell cycle transitions are 
present in comparable concentrations, we cannot use Michae- 
lis-Menten rate laws but must write out the mass-action rate 
laws for all species in the mechanism. On the other hand, 
our model is a little simpler than that of Markevich et al. 
[30], because we assume that the de-modifying enzyme E 
(comparable to their phosphatase) operates as a mass-action 
catalyst (E + IM — > I + E). Our mechanism is bistable even if 
E de-modifies IM and IMM equally rapidly. Because E is a 
mass-action catalyst, we can absorb its concentration into the 
rate constants, /c dml and k dm2 , for the de-modification reactions. 

3.2. SIMM*: a simplified motif for bistability at cell 
cycle transitions 

For the rate constant values given in table 2a, the Michaelis 
constant for the second modification, K m2 = fci2 + ^cat2)/ 
Ksi = 1-02 U, is large compared with K ml (0.01 U) for the 
first modification. (In our nomenclature, 'i<C m ' refers to a clas- 
sical Michaelis constant for an enzyme-catalysed reaction, 
and 'U' represents one unit of concentration of substrate.) 
Hence, it is reasonable to propose that the second modifi- 
cation can be modelled by a mass-action rate law in place 
of the Michaelis-Menten mechanism used so far. The simpli- 
fied motif (SIMM*) in figure 3a is described by a smaller set 



of ODEs: 



dJ 

dt ~ 
dt 

dJMM 
dt 

and 



-k as \A • I + kdilC + fcdml^Mj 
^catlC — k m A • I M — /Cdml^M + ^dm2^MM 5 



: k m A ■ I M — k dm2 h 



dC 
~dt 



dA 
"dt 



: k as \A • I — kdnC — /c ca tiC 



(3.1') 
(3.2') 
(3.3') 
(3.4') 



where the new rate constant is given by k m = k cat2 / K m2 = k cat2 
k as i/ (k ca t2 + k di2 ) ~ k as2 . (In our nomenclature, '/c m ' is a 
second-order rate constant for mass-action modification of 
substrate IM by enzyme A.) From the control plane for 
SIMM* (figure 3b), we see that the simplified motif repro- 
duces almost exactly the bistability properties of the full 
motif (figure 2d). 

Equations (3.1') -(3.4') satisfy two conservation conditions: 
A T = A + C = const, and J T = I + C + I M + J MM = const. 
As mentioned previously, because A T and J T are compara- 
ble, we cannot make a quasi-steady-state assumption on 
the enzyme -substrate complex, C, as in the derivation of the 
Michaelis-Menten rate law. However, as shown by Borghans 
et ah [31] and Ciliberto et ah [32], we may, under these circum- 
stances, use a 'total' quasi-steady-state assumption (TQSSA) to 
simplify equations (3.1') -(3.4') further. Following the TQSSA, 
we introduce a new variable, 1 = 1 + 0, and then solve for 
the steady state concentration of the A:I complex in terms of 
A T and I: K ml C = A • I = (A T - C)(J - C) = C 2 - (A T + J)C+ 
A T I, where K ml = (k dil + k catl )/k asl . 

Hence, 



C 



2A T I 



A T + I + X ml + ^/ (A T + J + X m i ) 2 - 4A T J 



(3.7) 



Because A T and X m i are constants, C is a function of the 
variable I. Now, equations (3.1 / )-(3.4 / ) reduce to a pair of 
nonlinear ODEs for the 'slow 7 variables I and I MM : 



dJ 

df : 

and 



= ^dml (It 
dlMM 

df 



I-lMM)-fccatlC(J) (3.8) 

k m (A T - C(f)) • (J T -I-Imm) - ^Imm- (3.9) 



In figure 3c, d we plot the nullclines of equation (3.8) and 
(3.9) in the (I, Imm) phase plane for two different values of 
J T - Note that these phase planes are identical to the sche- 
matic drawings in figure lc,d, because 1 = 1 + 0 is a 
measure of 'functional inhibitor' (the two forms of the 
inhibitor that are functional in neutralizing activator A) 
and I M m is a measure of 'functional activator' (the inhibitor 
is doubly modified only when the activator is functioning 
efficiently). 



3.3. Applying SIMM* to the phosphorylation of cyclin- 
dependent kinase inhibitors by cyclin-dependent 
kinases 

CKIs are an important class of CDK regulators in the cell 
cycle. CKIs bind to cyclin A- and cyclin B-dependent kinases 
to form inactive trimers (CKI:cyclin:Cdk). Because cyclin 
A:Cdk2 is a crucial initiator of DNA synthesis in mammalian 
cells, CKIs are important stabilizers of the 'late' Gl phase 



of the cell cycle (after the restriction point and before cells 
initiate a new round of DNA synthesis). The foremost 
example of a CKI is Sicl, a protein encoded by the SIC1 
gene in budding yeast [33]. The root name 'SIC stands 
for 'substrate and inhibitor of Cdk'. Sicl binds to and inhi- 
bits Clbx:Cdkl dimers [34], where the six B-type cyclins in 
budding yeast (Clbl-6) play essential roles in S/G2/M pro- 
gression. Sicl is also a target for multi-site phosphorylation 
by Clbx:Cdkl [22]. Hence, the interactions between Sicl- 
and Clb-dependent kinase have the topology of a SIMM 
motif. Sicl can also be phosphorylated by Cln-dependent 
kinases in budding yeast [17], a fact that will become 
important later. 

We base our new model (figure 4a) on the SIMM* motif in 
figure 3a, supplemented with reactions for the synthesis and 
degradation of Sicl. Because figure 4 refers specifically to 
budding yeast proteins, we change our notation somewhat: 
'CKI' becomes 'Sic' and 'cyclin:Cdk' dimers are referred to 
by their cyclin subunits 'Clb' and 'On'. With this notation, 
the mechanism in figure 4a is described by a further set 
of ODEs: 



d[SicT] 

dt 
d[SicP] 

dt 



■k s -k d [SicT] - k p _ h [Clb] [SicP] , (3.10) 

: ^catJb 

[cibSk] - (*a p + k d + k p _ h [Clb]) [SkP] 

(3.11) 



i d[Clb:SlC] = k as _ h [Clb] [Sic] - (k dLh + k d + k cath ) [Clb:Sic] , 



dt 



(3.12) 



where [Sic] = [SicT] - [Clb:Sic] - [SicP] and [Clb] = 
[ClbT] — [Clb:Sic]. The new rate constants are k s , the rate of 
synthesis of Sicl, and k d/ the 'background' rate of degradation 
of all forms of Sicl. We assume that the doubly phosphory- 
lated form of Sicl is rapidly degraded and thus lost 
immediately from the total pool of Sicl subunits. The mean- 
ings and assigned values of the rate constants in equations 
(3.10)-(3.12) are given in table 2c. 

Equations (3.10) -(3.12) comprise a bistable switch, as 
demonstrated in figure 4b. For low levels of total Clb-dependent 
kinase, the checkpoint is engaged (total Sicl level is high). As 
Clb accumulates in the cell, [ClbT] increases and, when it 
passes the threshold, Sicl is phosphorylated a second time 
and rapidly degraded (i.e. the checkpoint disengages). 

3.4. Extending the cyclin-dependent kinase -cyclin- 
dependent kinase inhibitor motif to the G1/S 
transition in budding yeast 

The model in the previous section provides a realistic mech- 
anism for bistability of the Gl/S transition in budding yeast. 
The two stable steady states of the bistable switch correspond 
to checkpoint engaged (Sicl abundant, Clb-kinase inactive) 
before the transition and to checkpoint disengaged (Sicl 
destroyed, Clb-kinase active) after the transition. In principle, 
the switch can be flipped by a signal that increases the total 
amount of Clb or by a signal that decreases the total 
amount of Sicl. In budding yeast, the signal comes from 
Cln-dependent kinase activity, which both upregulates syn- 
thesis of Clb proteins (by activating a transcription factor 
for Clb5 and Clb6) and downregulates Sicl (by multi-phos- 
phorylation and subsequent degradation). The Cln-type 
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Figure 4. The Sicl-cyclin switch in budding yeast, (a) Wiring diagram for Sid and Clb (B-type cyclins). Sid is a substrate and inhibitor of Clb-dependent kinase. 
The second phosphorylation of Sid leads to its rapid degradation, (b) Signal -response curve. The total concentration of Sid in the cell (a variable) depends, in a 
bistable manner, on the total concentration of Clb cyclins in the cell (a parameter, in this model). The black dashed line indicates how the switch flips from the 
Clb-inactive state to the Clb-active state as [ClbT] increases, (c) Wiring diagram for Sid, Clb and Cln (Gl-type cyclins). Sid is a substrate — but not an effective 
inhibitor — of Cln-dependent kinase, (d) Control plane, [ClnT] versus [ClbT]. The V-shaped bistable zone is coloured red. The 'cell cycle trajectory' (path a-b-c- 
d-e) is explained in the text, (e) Simulation of the Glb/S transition in budding yeast, and resetting the checkpoint at cell division. Equations (3.13) — (3.16) 
are solved numerically for the parameter values in table 2d, and with additional differential equations for the synthesis and degradation of ClnT and ClbT. 
The synthesis terms for ClnT and ClbT are turned on at t = 0 (the Start transition); at t = 30 min, synthesis of ClnT is turned off and degradation is turned 
on; at t = 50 min (cell division), synthesis of ClbT is turned off and degradation is turned up. The G1/S transition takes place at t w 20 min, when the cell 
cycle trajectory passes point b in (d), and there is enough active Clb-kinase to trigger DNA replication. About 50 min after Start, the cell divides and the G1 
phase is re-established. 



cyclins (Clnl and Cln2) are synthesized in late Gl. ClmCdkl 
dimers are not inhibited by Sicl, but they are capable of phos- 
phorylating Sicl. In this section, we couple the interactions of 
Sicl with Cln- and Clb-dependent kinases to build a model of 
the Gl/S transition in budding yeast. 

The model (figure 4c) is composed of two modules: a 
Clb module (from the previous section), where Sicl binds 
strongly to Clb:Cdkl dimers to form inactive trimers (i.e. 
enzyme -substrate complexes that only slowly convert to 
enzyme + SiclP); and a Cln module, where Sicl is rapidly 



phosphorylated by Cln:Cdkl. The model in figure 4c is 
described by the following ODEs: 



d[SicT] 
dt 



: k s - /c d [SicT] - (*pjb[Clb] + fc p ^[Cln])[SicP], 



d [Clb: Sic] 
dT~ 



= k ash [Sic] [Clb] - (k dLh + k d + k cath ) [Clb:Sic] , 



(3.13) 



(3.14) 



d[Cln:Sic] 



dt 



-- k asjl [Sic] [Cln] - (k disi + k d + k caisl ) [ClnSic] 



(3.15) 



d[SicP] 



and ^ = k cat _ h [ClbSic] + /c cat _n [Cln:Sic] - (k dp + k d ) 
x [SicP] - (fc p _ b [Clb] +/c p _ n [Cln])[SicP], 

(3.16) 

where [Sic] = [SicT] - [Clb:Sic] - [Cln:Sic] - [SicP], [Clb] = 
[ClbT] - [Clb:Sic] and [Cln] = [ClnT] - [Cln:Sic]. All calcu- 
lations for this model are done for the rate constant values 
given in table 2d. 

Figure 4d shows the region of bistability in the control 
plane spanned by [ClnT] and [ClbT]. To explain irreversibil- 
ity of the Gl/S transition, we plot a schematic 'cell cycle 
trajectory 7 (path a-b-c-d-e) on the control plane. During 
the second half of Gl phase (from Start to the onset of S 
phase), both [ClbT] and [ClnT] are increasing, as indicated 
by the dashed line a-c. At the beginning of this process 
(from a to b), the checkpoint mechanism is engaged (Sicl 
abundant and Clb-kinase inactive). As the trajectory passes 
point b, the checkpoint mechanism disengages (Sicl is 
destroyed and Clb-kinase is abundant and active). After the 
Gl/S transition, high activity of Clb-kinase shuts off the pro- 
duction of Clnl and Cln2, and [ClnT] drops along path c, d; 
however, along this path, [ClbT] is large enough to keep the 
checkpoint mechanism disengaged. Only much later in the 
cell cycle (during mitosis) are the Clb cyclins destroyed 
(path d-e-a), and the checkpoint re-engages (at point e). 
In figure 4e we simulate this scenario, to show that our 
model of the Gl/S transition does indeed disengage and 
re-engage as described. 

Our model confirms that both the bistability of the Gl/S 
transition and the abruptness of Sicl degradation are depen- 
dent on Clb cyclins [22]. These conclusions are consistent 
with kinetic data on Sicl degradation in individual budding 
yeast cells (C. Tang 2012, personal communication). 

3.5. Applying SIMM* to the mitotic checkpoint 

The mitotic checkpoint (also known as the spindle assembly 
checkpoint) guards the metaphase-to-anaphase transition. 
As cells enter mitosis, each chromosome is fully replicated, 
and sister chromatids are held together by cohesin rings. 
Each duplicated chromosome has a centromeric region with 
kinetochores on each sister chromatid. During prometaphase, 
every chromosome must become attached by its sister kineto- 
chores to opposite poles of the mitotic spindle. The mitotic 
checkpoint delays progression into anaphase (sister chroma- 
tid separation) until biorientation of the chromosomes is 
complete. When the last chromosome pair comes into align- 
ment on the spindle, the checkpoint is disengaged, and 
separation of sister chromatids is initiated by cleavage of 
cohesin rings, catalysed by separase [35]. Separase is kept 
inactive by its stoichiometric binding partner, securin. In 
addition, high activity of cyclin B:Cdkl in prometaphase 



and metaphase inhibits further mitotic progression [36]. 
Both securin and cyclin B are removed by APC-dependent 
poly-ubiquitination followed by proteasome-catalysed degra- 
dation [37,38]. APC requires the Cdc20 cof actor in order to 
ubiquitinate securin and cyclin B during mitotic progression. 
Therefore, APC:Cdc20 is an activator of mitotic progression 
that drives the cell into anaphase. 

The mitotic checkpoint has three essential dynamical fea- 
tures: (i) even one unattached kinetochore is sufficient to keep 
the checkpoint engaged [39]; (ii) after the last kinetochore 
becomes attached, the checkpoint rapidly disengages [40]; 
and (iii) in anaphase, when the chromosome alignment 
signal drops to zero, the checkpoint does not re-engage 
[41-43]. These features suggest that dynamical changes in 
the checkpoint mechanism are crucial to proper progression 
through metaphase and anaphase. We have proposed that 
the essential features of the mitotic checkpoint derive from 
its bistable switching properties [44]. 

Unattached or mis-attached kinetochores catalyse the for- 
mation of the MCC, which binds strongly and reversibly to 
the APC that inhibits APC binding to free Cdc20. The 
APC:MCC complex includes many proteins, including 
Mad2, Bub3, Mad3/BubRl and Cdc20, as well as the APC 
core complex [45]. As required by the SIMM motif, MCC is 
a substrate as well as an inhibitor of APC. MCC inactivation 
is an active, energy-requiring process catalysed by APC [23]. 
The Cdc20 subunit of the APC:MCC complex is poly-ubiqui- 
tinated by APC, which causes degradation of Cdc20 and 
disassembly of the complex [23,25-27,46,47]. Because 
Cdc20 level is not limiting during M phase, new MCCs are 
assembled as long as unattached kinetochores are present, 
giving a signal that prometaphase is still in progress. When 
a chromosome becomes bioriented, its centromeric region 
comes under tension (being pulled in opposite directions by 
the kinetochore attachments), and this tension silences kine- 
tochore signalling to the MCC. 

In summary, MCC and APC are an inhibitor and an 
activator locked in a double-negative feedback loop [48], 
as required for a SIMM motif (figure 5a). From the reac- 
tion mechanism in figure 5a, it is obvious that the total 
concentration of APC, [APCT] = [APC] + [APC:MCC], is 
constant and that the total concentration of Mad2, 
[Mad2T] = [Mad2] + [MCCU] + [APC:MCC] + [MCC], is 
also constant. It is convenient to define a new variable 
[MCCT] = [MCC] + [APC:MCC] + [MCCU] = total concen- 
tration of mitotic checkpoint complexes. Using these 
definitions, we can write a pair of nonlinear ODEs that describe 
our model of the interactions between MCC and APC: 



d[MCCT] 
dt 



: fc a N T (l - X tens )([M2T] - [MCCT]) 
- k u [APC] [MCCU] , 



(3.17) 



d[MCCU] = ^ t [APC:MCC] - (k du +k u [APC]) [MCCU], 



dt 



(3.18) 



a? 

o' 



and 



[APC] = 



[APCT] + [MCCU] - [MCCT] - K n 



([APCT] + [MCCU] - [MCCT] - K m ) 2 + 4K m [APCT] 



(3.19) 
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Figure 5. The mitotic checkpoint, (a) Wiring diagram. The mitotic checkpoint complex (MCC) binds to and inhibits the anaphase-promoting complex (APC). 
APC ubiquitinates MCC, and the ubiquitin moiety is removed by a de-ubiquitinase. Poly-ubiquitination of MCC leads to degradation of some of its components 
(notably Cdc20) and release of inactive checkpoint proteins (labelled by Mad2 only). Checkpoint proteins (Mad2, etc.), are reactivated by tensionless centromeres, 
N 0 = /V T (1 - /tens), leading to reassembly of MCC. Securin degradation by active APC leads to dissolution of cohesin complexes and separation of sister chromatids 
in early anaphase, (b) Signal -response curve: [MCCT], the total concentration of active MCC, versus X tensi the fraction of centromeric regions that are in tension on 
the mitotic spindle. The black dashed curve indicates that the mitotic checkpoint is engaged (i.e. [MCCT] large) until the last chromosome aligns on the mitotic 
spindle. When / tens increases above approximately 0.97, the checkpoint disengages and APC is activated. As the cell enters anaphase, X Xem drops back to zero, (c) 
Simulation of in vitro release of the mitotic checkpoint. Compare with fig. 1 of Reddy et al. [23]. Note that when X lens drops below approximately 0.89 the check- 
point re-engages, as predicted by the signal-response curve in (b). (d) Wiring diagram of the CycB-MCC-APC network. Cycin B-dependent kinase is required for 
activation of tensionless centromeres, X 0k . (e) Signal -response curve: [MCCT] versus X tem . Now the threshold for re-engaging the checkpoint has moved to negative 
values of X temi and the checkpoint remains disengaged when / tens drops to zero in anaphase (black dashed line), (f) Simulation of in vivo release of the mitotic 
checkpoint. The checkpoint mechanism does not re-engage as / tens decreases. 
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where [APC:MCC] = [APCT] - [APC] and K m = (k di + 
k cat )/k as . Equation (3.19) is derived, according to the TQSSA, 
by solving the following algebraic equation for [APC]: 
(^di+^cat)([APCT]-[APC])-^ as [APC]([MCCT] - [MCCU]- 
[APCT] + [APC]). 

An equation for APC-mediated degradation of securin 
is added as an indicator of the metaphase-to-anaphase 
transition: 

drSecl 

= k ssec - (k dsec + fc d sec,apc[APC])[Sec]. (3.20) 

The parameters N T and X tens in equation (3.17) are the 
total number of chromosomes in the cell and the fraction of 
centromeres that are under tension because their two kineto- 
chores are attached to opposite poles of the mitotic spindle. 
We assume that Cdc20 is present in excess during mitosis 
and MCC formation is limited by Mad2 (or any other MCC 
subunit). However, as long as APC is occupied by MCC it 
cannot ubiquitinate its substrates (Securin and cyclin B), 
whose degradation is required for anaphase progression. 
During prometaphase, X tens increases from 0 to 1. The mean- 
ings and assigned values of the other constants in equations 
(3. 17) -(3.20) are given in table 2e. 

At the beginning of prometaphase, when X tens = 0, the 
control system is trapped in a stable steady state with lots 
of active MCC, and most of the APC is bound up in inactive 
APC:MCC complexes. As chromosomes come into alignment 



on the spindle, the checkpoint remains engaged until X tens 
gets very close to 1 (figure 5b). For the parameter set in 
table 2e, the mitotic checkpoint does not disengage until 

Xtens ^ > Xdisengage ~ 0.97. 

The evident 'sensitivity' of the mitotic checkpoint to even 
a single unaligned chromosome is not a freak consequence of 
our choice of rate constants in table 2e. In the electronic sup- 
plementary material, §S2 we use two-parameter bifurcation 
diagrams to show that 0.95 < X disengage < 1 over a broad 
range of values of the rate constant in equations (3. 17) -(3. 19). 

In figure 5c we show a simulation of how fast the MCC dis- 
engages after the last centromere comes under tension. At first 
there is a considerable delay, as the excess supply of 'func- 
tional' MCC, [MCC] + [APC:MCC], is slowly ubiquitinated 
by the APC. After approximately 45 min, in our simulation, 
total functional MCC becomes comparable with total APC, 
and there is now a large enough concentration of free APC to 
carry out the second stage of ubiquitination, which is an 
accelerating, autocatalytic process in our simulations and 
in experimental measurements [23]. In vitro experiments on 
the processing of APC MCC complexes show a long delay 
(approx. 30 min) before anaphase substrates (securin and 
cyclin B) start to be destroyed [23], as predicted by our 
model. However, the long delay in figure 5c is not consistent 
with in vivo observations [40], where anaphase substrates are 
degraded within minutes of lifting the mitotic checkpoint. 
This disagreement between model and experiment suggests 



that some other processes are at work in vivo. The long delay in 
figure 5c is a result of the rapid rate of de-ubiquitination of 
MCCU, which is a necessary requirement for sensitivity of 
the mitotic checkpoint. During the lag phase, as functional 
MCC is ubiquitinated by APC, the de-ubiquitinating enzyme 
(USP44 in mammalian cells) is rapidly recycling MCCU to 
MCC. Hence, it takes a long time to poly-ubiquitinate MCC 
and destroy the inhibitor. Seeing it this way, we suggest that 
in vivo there are mechanisms that downregulate the activity 
of the de-ubiquitinase as the mitotic checkpoint disengages, 
so that functional MCC can be more rapidly destroyed. 

3.6. Extending the anaphase-promoting complex - 
mitotic checkpoint complex motif to account 
for irreversibility of the metaphase-to- 
anaphase transition 

In §3.5 we showed how a SIMM* motif might account for 
'sensitivity' of the mitotic checkpoint and how it might be 
modified to account for 'fast release'. In this section, we 
describe an extension of the model to account for 'irreversibil- 
ity' of the metaphase-to-anaphase transition. 

The signal -response diagram in figure 5b shows that the 
mitotic checkpoint disengages when the fraction of bioriented 
chromosomes gets close to 1. After the transition, centromeric 
cohesins are destroyed and sister chromatids move towards 
opposite poles of the spindle. Because their kinetochores 
are no longer silenced by centromeric tension, we might 
expect the mitotic checkpoint signal to go back to a large 
value. If so, then figure 5b predicts that the checkpoint 
would re-engage when X tens drops below approximately 
0.89. To account for irreversibility of the metaphase-to- 
anaphase transition, we must introduce additional regulatory 
steps to prevent re-activation of the MCC during anaphase. 

Reactivation of the mitotic checkpoint in anaphase does 
not happen during normal mitotic progression because 
cyclin B:Cdkl activity appears to be essential for checkpoint 
signalling [49]. Cyclin B:Cdkl activity is high in prometa- 
phase when tensionless kinetochores are signalling. But 
cyclin B is degraded by the APC during metaphase, so the 
mitotic checkpoint cannot be reactivated during anaphase. 
The non-resettable character of the mitotic checkpoint 
depends on cyclin B degradation by APC [42,43]. Therefore, 
mitotic checkpoint control is based on one activator (APC) 
and two inhibitors (cyclin B:Cdkl and MCC). In addition to 
the SIMM motif involving APC and MCC, APC and cyclin 
B are engaged in a double-negative feedback loop: (i) cyclin 
B:Cdkl indirectly inactivates APC by promoting MCC 
assembly at tensionless kinetochores; and (ii) cyclin B 
becomes a substrate of APC-mediated ubiquitination and 
degradation, once the checkpoint is silenced (figure 5d). 

The extended SIMM motif (figure 5d) is governed by the 
following ODEs: 

d[C d y f CB] = k scyc - (k dcyc + /c dcyc _ apc [APC])[CycB], (3.21) 



d[MCCU] 
dt 



£ cat [APC:MCC] - (k du + fc u [APC]) [MCCU] 



Xqa = MCycBKl - X tens - X 0A ) - £ in [CAPP]X 0A , 



dt 

d[MCCT] 
dt 



(3.22) 



and d[AP ^; MCC1 = k as [APC] [MCC] - (k di +k cat ) [APC : MCC] 



(3.24) 
vlCC], 
(3.25) 



: fc a N T Xo A ([M2T] - [MCCT]) - k u [APC] [MCCU] , 

(3.23) 



where [APC] = [APCT] - [APC:MCC] and [MCC] = 
[MCCT] - [APC:MCC] - [MCCU]; and X 0A indicates the 
fraction of centromeres without tension that are activating 
the MCC. Given the parameter values in table If, this 
model predicts that for X tens = 0 the mitotic checkpoint has 
two stable steady states (figure 5e): 'engaged' ([MCCT] 
high, [CycB] high, [APC] low) and 'disengaged' ([MCCT] 
low, [CycB] low, [APC] high). The checkpoint starts in the 
engaged state in prometaphase with X tens = 0 and ends up 
in the disengaged state in telophase with X tens = 0. 

One might well ask, if the mitotic checkpoint is inactive 
in the final state of the one-way switch (figure 5e), how 
does it get activated again in the next cell cycle? During 
Gl phase of the next cycle, the mitotic activity of APC 
(which is dependent on an auxiliary protein, Cdc20) is lost, 
because the Cdkl /CycB-dependent production of Cdc20 is 
terminated, Cdc20 protein is rapidly degraded and Cdkl/ 
CycB-dependent activatory phosphorylations of APC are 
lost. The inactivation of APC/Cdc20 allows reactivation of 
the mitotic checkpoint in Gl phase of the next cycle. 

In the electronic supplementary material, §S3 we investi- 
gate the robustness of this model of the mitotic checkpoint 
with respect to the dual requirements of sensitivity and irre- 
versibility. We show that the region of parameter space where 
equations (3.21) -(3.25) exhibit both properties is quite 
restricted, compared with the robustness of our earlier 
model of the mitotic checkpoint [44]. The major difference 
between the present model and that of He et al. [44] is that 
the counter-acting protein phosphatase (CAPP in equation 
(3.22)) is regulated by CycB-kinase in [44], but not here. 
Hence, we may conclude that the simple model considered 
here is sufficient to account for sensitivity and irreversibility 
of the mitotic checkpoint. However, the mechanism is not 
very robust and release from the checkpoint is rather sluggish. 
The robustness issue can be solved by additional regulatory 
signals on the counter-acting protein phosphatase, and the 
sluggishness issue can be solved by additional regulatory 
signals on the counter-acting de-ubiquitinase enzyme. 



4. Discussion 

4.1. An underlying design principle for eukaryotic cell 
cycle transitions 

Cell cycle transitions are impelled by specific activators: Cdks 
drive the transitions from Gl to S and from G2 to M, whereas 
the APC drives exit from mitosis back to Gl. In the pre-tran- 
sition state, the activator is shut down by an inhibitor. The 
inhibitors and activators functioning at specific cell cycle tran- 
sitions are listed in table 1. In each case, prior to the transition, 
the inhibitor binds to the activator and inactivates it, either 
stoichiometrically (by forming an inactive heteromeric com- 
plex) or catalytically (by covalently modifying the activator). 
After the transition, the activator inactivates the inhibitor, 



either directly or indirectly. This double-negative feedback 
loop between activator and inhibitor is the basis of a bistable 
switch that makes the cell cycle transition a unidirectional 
(irreversible) change. 

In this work, we analyse the case where the activator 
inactivates the inhibitor by direct modification. In this case, 
the inhibitor is both a stoichiometric inhibitor and a sub- 
strate of the activator. This relationship between inhibitor 
and activator is exhibited in the cell cycle at the Gl/S transition 
(activator = Cdk and inhibitor = CKI) and at the meta-to- 
anaphase transition (activator = APC, inhibitor = MCC). We 
show that the antagonism between inhibitor and activator cre- 
ates a bistable switch if (i) the inhibitor is modified by the 
activator more than once (multi-site modification) in separate 
collisions (distributive mechanism) and (ii) the modifications 
of the inhibitor are reversed by counter-acting reactions. 

Both of these requirements are satisfied at the Gl/S and 
meta/ anaphase transitions. CKIs get phosphorylated at mul- 
tiple sites, while the Cdc20 component of the MCC gets poly- 
ubiquitinated. CKI phosphorylation and MCC ubiquitination 
are reversed by phosphatase and de-ubiquitinase enzymes, 
respectively. The bistable switch is characterized by two 
sharp thresholds. The thresholds are sensitively dependent 
on upstream regulators and the counter-acting enzymes 
(the phosphatase or de-ubiquitinase) in the network. 
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Figure 6. Wiring diagrams for other cell cycle transitions, (a) The G2/M tran- 
sition. Wee1 is a protein kinase that phosphorylates and inactivates Cdk1 (i). 
Simultaneously, active Cdkl (complexed with cyclin B) phosphorylates Wee1 
on multiple sites in a typical SIMM motif (ii), where Wee1 functions as a 
stoichiometric inhibitor of Cdkl as well as a substrate. Cdc25, the phosphatase 
that dephosphorylates P-Cdkl, is activated by multi-site phosphorylation cat- 
alysed by active Cdkl. (b) Mitotic exit. Cdhl (in combination with the APC) 
stabilizes G1 phase of the cell cycle by poly-ubiquitinating cyclin B (i). Sim- 
ultaneously, active CdkkCycB phosphorylates Cdhl on multiple sites (ii). Both 
sides of the interaction have SIMM topology, but it is not known if either of 
the 'substrates' serves also as a 'stoichiometric inhibitor' of the enzyme. 
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4.2. Application of this principle to other cell cycle 
transitions and checkpoint mechanisms 

Our proposed network can be easily extended to the Start 
transition in yeasts (and the analogous restriction point in 
mammalian cells). The activator of the Start transition is a 
transcription factor, SBF, that drives expression of Gl cyclins, 
Clnl and Cln2. (In mammalian cells, the analogous com- 
ponents are the E2F family of transcription factors and the 
family of E-type cyclins.) The inhibitor is Whi5 (the RB in 
mammalian cells), which binds to and inactivates SBF. 
In this case, inactivation of the inhibitor by the activator is 
indirect: the SBF-induced, Gl cyclin-dependent kinases phos- 
phorylate and inactivate the inhibitor, Whi5. It is noteworthy 
that Whi5 is multiply phosphorylated by CdkliCln com- 
plexes. Hence, the mechanism of the Start transition bears 
many similarities to our SIMM motif. 

The G2/M transition in eukaryotes is promoted by 
CdkliCycB complexes, which are kept inactive in G2 phase 
by inhibitory phosphorylation of Cdkl subunits, catalysed 
by Weel-like kinases (figure 6a). Weel (the inhibitor of our 
scheme) gets phosphorylated at multiple sites by CdkliCycB 
(the activator), and multiply phosphorylated Weel is less 
active [12]. This case differs from our SIMM motif in that 
the inhibitor inactivates the activator by covalent modifi- 
cation after complex formation. But stoichiometric inhibition 
of CdkliCycB by Weel is also at play here, as indicated by 
the phenotype of fission yeast cells carrying mutant Cdkl 
subunits that cannot be phosphorylated by Weel. These 
mutant cells (cdkl AF ) are poorly viable, whereas fission 
yeast cells lacking both inhibitory kinases (weel ts miklA) are 
totally inviable. The double mutant cells undergo mitotic cat- 
astrophe (i.e. they enter mitosis before they complete DNA 
replication) and die, whereas some mutant cells carrying 
non-phosphorylable Cdkl subunits survive, presumably 
because the inhibitory kinases delay the G2/M transition 



solely on the basis of stoichiometric (rather than catalytic) 
inhibition of CdkliCycB. The advantage of enzymatic inhi- 
bition of the activator at the G2/M transition is that much 
less inhibitor (Weel) is required to stabilize the G2 state. 
The disadvantage is that the mechanism requires a second 
enzyme (Cdc25 phosphatase) to reverse this modification. 
Interestingly, Cdc25 is also phosphorylated by CdkliCycB 
at multiple sites; these phosphorylations are activatory, creat- 
ing a positive feedback loop between Cdkl and Cdc25. 

Finally, in Gl phase of the budding yeast cell cycle, the 
APC combines with a different auxiliary protein, Cdhl, to 
serve as an inhibitor of the Gl/S transition and an activator 
of the mitotic exit transition. APOCdhl interacts with 
CdkliCycB (an activator of the Gl/S transition, and an 
inhibitor of mitotic exit). APCiCdhl inhibits CdkliCycB by 
poly-ubiquitinating CycB, which induces CycB degradation 
by the proteasome. CdkliCycB inhibits APC /Cdhl by phos- 
phorylating the Cdhl subunit on multiple sites. Hence, this 
regulatory motif consists of two SIMMs back to back 
(figure 6b). In this topology, bistability does not require 
that either substrate be a stoichiometric inhibitor of the 
enzyme to which it binds. However, bistability does require 
multiplicity of one or both of the catalytic processes 
(poly-ubiquitination and /or multiple phosphorylation). 

As simple, robust mechanisms for bistable, irreversible 
transitions, SIMM motifs may underlie other cellular 
'decision-making' processes in signalling networks and 
developmental pathways. 

5. Material and methods 

5.1. Mathematical modelling and simulation 

Each cell cycle transition mechanism is represented as a 
wiring diagram, which is then converted into a set of rate 



equations (nonlinear ODEs) by the law of mass action. The 
ODEs are solved numerically using a software package 
called XPPAUT, which is freely available from http:// 
www.math.pitt.edu/~bard/xpp/xpp.html. Before doing 
any simulations, we must assign numerical values to all 
the rate constants introduced by the law of mass action. 
We make these assignments intuitively, according to the 
following principles: 

— The rate constants should be 'reasonable' numbers. For 
example, the rate of association of two proteins to form 
a complex should not be faster than the diffusion-limited 
rate for second-order chemical reactions, which (for pro- 
teins in cytoplasm) is approximately 10 9 1 mol -1 s -1 , or 
approximately 600 U -1 min -1 in our units. The degra- 
dation rates of proteins such as cyclin and CKI should 
be consistent with observed rates (half-life~min), and 
their rates of synthesis should be consistent with intra- 
cellular concentrations in the 10-100 nM range. 

— The rate constant values should give signal -response 
curves (one-parameter bifurcation diagrams) and control 
planes (two-parameter bifurcation diagrams) that are con- 
sistent with the observed dynamic properties of cell cycle 
checkpoints (bistability, irreversibility and fast release). 



In addition to providing algorithms for numerical solutions 
of ODEs, XPPAUT provides a suite of tools for the analysis of 
dynamical systems, including phase plane portraits and bifur- 
cation analysis. We use these tools to illustrate the dynamical 
properties of checkpoint mechanisms. For example, figure 3c 
and 3d are phase plane portraits of the SIMM* motif, the 
'signal -response curves' in figures 2, 4 and 5 are one-parameter 
bifurcation diagrams, and the 'control planes' in figures 2, 3 and 
4 are two-parameter bifurcation diagrams. Interested readers 
should consult the electronic supplementary material, §S4 for 
a simple introduction to these tools. 

All figures in this paper can be reproduced in XPPAUT 
using the '.ode' files provided in the electronic supplemen- 
tary material, §S5. 

6. Acknowledgements 

This work was supported in part by the National Science Foun- 
dation's award number DBI-0904340 (to A.V.), by the National 
Institutes of Health's award no. 5R01GM078989-06 (to J.J.T.), 
and by the BBSRC and European Community's Seventh 
Framework Programmes UniCellSys/201142 and MitoSys/ 
241548 (to B.N.). 



References 

1. Nasmyth K, Hunt T. 1993 Dams and sluices. Nature 
366, 634-635. (doi:10.1038/366634a0) 

2. Nasmyth K. 1995 Evolution of the cell cycle. Phil. 
Trans. R. Soc. Lond. B 349, 271-281. (doi:10.1098/ 
rstb.1995.0113) 

3. Novak B, Tyson JJ. 1993 Numerical analysis of a 
comprehensive model of M-phase control in 
Xenopus oocyte extracts and intact embryos. J. Cell 
Sci. 106, 1153-1168. 

4. Tyson JJ, Novak B, Chen K, Val J. 1995 Checkpoints in 
the cell cycle from a modeler's perspective. Prog. Cell 
Cycle Res. 1, 1-8. (doi: 1 0.1 007/978-1-461 5-1 809-9_1 ) 

5. Novak B, Tyson JJ, Gyorffy B, Csikasz-Nagy A. 2007 
Irreversible cell-cycle transitions are due to systems- 
level feedback. Nat. Cell Biol. 9, 724-728. (doi:10. 
1038/ncb0707-724) 

6. Cross FR, Archambault V, Miller M, Klovstad M. 2002 
Testing a mathematical model of the yeast cell 
cycle. Mol. Biol. Cell 13, 52-70. (doi:10.1091/mbc. 
01-05-0265) 

7. Lopez-Aviles S, Kapuy 0, Novak B, Uhlmann F. 2009 
Irreversibility of mitotic exit is the consequence of 
systems-level feedback. Nature 459, 592-595. 
(doi : 1 0. 1 038/nat u re07984) 

8. Pomerening JR, Sontag ED, Ferrell Jr JE. 2003 
Building a cell cycle oscillator: hysteresis and 
bistability in the activation of Cdc2. Nat. Cell Biol. 5, 
346-351. (doi:10.1038/ncb954) 

9. Sha W, Moore J, Chen K, Lassaletta AD, Yi CS, Tyson JJ, 
Sible JC. 2003 Hysteresis drives cell-cycle transitions in 
Xenopus laevis egg extracts. Proc. Natl Acad. Sci. USA 
100, 975-980. (doi:10.1073/pnas.0235349100) 

10. Yao G, Lee TJ, Mori S, Nevins JR, You L. 2008 A 
bistable Rb-E2F switch underlies the restriction point. 
Nat. Cell Biol. 10, 476 - 482. (doi:10.1038/ncb1711) 



11. Holt LJ, Krutchinsky AN, Morgan DO. 2008 Positive 
feedback sharpens the anaphase switch. Nature 
454, 353-357. (doi:10.1038/nature07050) 

12. Kim SY, Ferrell Jr JE. 2007 Substrate competition as a 
source of ultrasensitive in the inactivation of Wee1. 
Cell 128, 1133-1145. (doi:10.1016/j.cell.2007.01.039) 

13. Zetterberg A, Larsson 0, Wiman KG. 1995 What 
is the restriction point? Curr. Opin. Cell Biol. 7, 
835 - 842. (doi:1 0.1 01 6/0955-0674(95)80067-0) 

14. Costanzo M et al. 2004 CDK activity antagonizes Whi5, 
an inhibitor of G1/S transcription in yeast. Cell 117, 
899-913. (doi:10.1016/j.cell.2004.05.024) 

15. de Bruin RA, McDonald WH, Kalashnikova Tl, Yates 
lllrd J, Wittenberg C. 2004 Cln3 activates Gl-specific 
transcription via phosphorylation of the SBF bound 
repressor Whi5. Cell 117, 887-898. (doi:10.1 016/j. 
cell.2004.05.025) 

16. Morgan DO. 2007 The cell cycle: principles of control. 
London, UK: New Science Press. 

17. Verma R, Annan RS, Huddleston MJ, Carr SA, 
Reynard G, Deshaies RJ. 1997 Phosphorylation of 
Siclp by G1 Cdk required for its degradation and 
entry into S phase. Science 278, 455-460. (doi:10. 
1126/science.278.5337.455) 

18. Benito J, Martin-Castellanos C, Moreno S. 1998 
Regulation of the G1 phase of the cell cycle by 
periodic stabilization and degradation of the 
p25rum1 CDK inhibitor. EM BO J. 17, 482-497. 
(doi:10.1093/emboj717.2.482) 

19. Sheaff RJ, Groudine M, Gordon M, Roberts JM, 
Clurman BE. 1997 Cyclin E-CDK2 is a regulator of 
p27Kip1. Genes Dev. 11, 1464-1478. (doi:10.1 101/gad. 
11.11.1464) 

20. Nash P, Tang X, Orlicky S, Chen Q, Gertler FB, 
Mendenhall MD, Sicheri F, Pawson T, Tyers M. 2001 



Multisite phosphorylation of a CDK inhibitor sets a 
threshold for the onset of DNA replication. Nature 
414, 514-521. (doi : 1 0. 1 038/35 1 07009) 

21. Chen KC, Csikasz-Nagy A, Gyorffy B, Val J, Novak B, 
Tyson JJ. 2000 Kinetic analysis of a molecular model 
of the budding yeast cell cycle. Mol. Biol. Cell 11, 
369-391. 

22. Koivomagi M, Valk E, Venta R, lofik A, Lepiku M, 
Balog ER, Rubin SM, Morgan DO, Loog M. 2011 
Cascades of multisite phosphorylation control Sid 
destruction at the onset of S phase. Nature 480, 
128-131. (doi:10.1038/nature10560) 

23. Reddy SK, Rape M, Margansky WA, Kirschner MW. 
2007 Ubiquitination by the anaphase-promoting 
complex drives spindle checkpoint inactivation. 
Nature 446, 921-925. (doi:10.1038/nature05734) 

24. Nilsson J, Yekezare M, Minshull J, Pines J. 2008 The 
APC/C maintains the spindle assembly checkpoint 
by targeting Cdc20 for destruction. Nat. Cell Biol. 
10, 1411-1420. (doi:10.1038/ncb1799) 

25. Varetti G, Guida C, Santaguida S, Chiroli E, 
Musacchio A. 2011 Homeostatic control of mitotic 
arrest. Mol. Cell 44, 710-720. (doi:1 0.1 016/j. 
molcel.201 1.1 1.014) 

26. Foster SA, Morgan DO. 2012 The APC/C subunit 
Mnd2/Apc15 promotes Cdc20 autoubiquitination 
and spindle assembly checkpoint inactivation. 
Mol. Cell 47, 921-932. (doi:1 0.1 01 6/j.molcel. 
2012.07.031) 

27. Uzunova K et al. 2012 APC15 mediates CDC20 
autoubiquitylation by APC/C(MCC) and disassembly 
of the mitotic checkpoint complex. Nat. Struct. Mol. 
Biol. 19, 1116-1123. (doi:10.1038/nsmb.2412) 

28. Salazar C, Hofer T. 2009 Multisite protein 
phosphorylation: from molecular mechanisms to 



kinetic models. FEBS J. 276, 3177-3198. (doi:10. 
1111/j.1742-4658.2009.07027.x) 

29. Ortega F, Garces JL, Mas F, Kholodenko Bl\l, Cascante 
M. 2006 Bistability from double phosphorylation 

in signal transduction: kinetic and structural 
requirements. FEBS J. 273, 3915-3926. (doi:10. 
1 1 1 1/j.l 742-4658.2006.05394.X) 

30. Markevich Nl, Hoek JB, Kholodenko BN. 2004 
Signaling switches and bistability arising from 
multisite phosphorylation in protein kinase 
cascades. J. Cell Biol. 164, 353-359. (doi:10.1083/ 
jcb.200308060) 

31. Borghans JA, de Boer RJ, Segel LA. 1996 Extending 
the quasi-steady state approximation by changing 
variables. Bull. Math. Biol. 58, 43-63. (doi:10. 
1007/BF02458281) 

32. Ciliberto A, Capuani F, Tyson JJ. 2007 Modeling 
networks of coupled enzymatic reactions using 
the total quasi-steady state approximation. PLoS 
Comput. Biol. 3, e45. (doi:10.1371/journal. 
pcbi.0030045) 

33. Nugroho TT, Mendenhall MD. 1994 An inhibitor of 
yeast cyclin-dependent protein kinase plays an 
important role in ensuring the genomic integrity of 
daughter cells. Mol. Cell Biol. 14, 3320-3328. 

34. Schwob E, Bohm T, Mendenhall MD, Nasmyth K. 
1994 The B-type cyclin kinase inhibitor p40SIC1 
controls the G1 to S transition in 5. cerevisiae. 
Cell 79, 233-244. (doi:1 0.1 01 6/0092-8674(94) 
90193-7) 



35. Nasmyth K, Peters JM, Uhlmann F. 2000 Splitting 
the chromosome: cutting the ties that bind sister 
chromatids. Science 288, 1379-1385. (doi:10.1 1 26/ 
science.288.5470.1379) 

36. Stemmann 0, Zou H, Gerber SA, Gygi SP, Kirschner 
MW. 2001 Dual inhibition of sister chromatid 
separation at metaphase. Cell 107, 715-726. 
(doi:1 0.1 01 6/S0092-8674(01 )00603-1 ) 

37. Peters JM. 2006 The anaphase promoting complex/ 
cyclosome: a machine designed to destroy. Nat. 
Rev. Mol. Cell Biol. 7, 644-656. (doi:10.1038/ 
nrm1988) 

38. Pines J. 2011 Cubism and the cell cycle: the many 
faces of the APC/C. Nat. Rev. Mol. Cell Biol. 12, 
427-438. (doi:10.1038/nrm3132) 

39. Rieder CL, Schultz A, Cole R, Sluder G. 1994 
Anaphase onset in vertebrate somatic cells is 
controlled by a checkpoint that monitors sister 
kinetochore attachment to the spindle. J. Cell Biol. 
127, 1301-1310. (doi:10.1083/jcb.127.5.1301) 

40. Clute P, Pines J. 1999 Temporal and spatial control 
of cyclin B1 destruction in metaphase. Nat. Cell Biol. 
1, 82-87. (doi:10.1038/10049) 

41. Parry DH, Hickson GR, O'Farrell PH. 2003 Cyclin B 
destruction triggers changes in kinetochore behavior 
essential for successful anaphase. Curr. Biol. 13, 
647-653. (doi: 1 0.1 01 6/S0960-9822(03)00242-2) 

42. Mirchenko L, Uhlmann F. 2010 Sli15(INCENP) 
dephosphorylation prevents mitotic checkpoint 
reengagement due to loss of tension at anaphase 



onset. Curr. Biol. 20, 1396-1401. (doi:10.1016/j. 
cub.2010.06.023) 

43. Oliveira RA, Hamilton RS, Pauli A, Davis I, Nasmyth 
K. 2010 Cohesin cleavage and Cdk inhibition trigger 
formation of daughter nuclei. Nat. Cell Biol. 12, 
185-192. (doi:10.1038/ncb2018) 

44. He E, Kapuy 0, Oliveira RA, Uhlmann F, Tyson JJ, 
Novak B. 2011 System-level feedbacks make the 
anaphase switch irreversible. Proc. Natl Acad. Sci. 
USA 108, 10 016-10 021. (doi:10.1073/pnas. 
1102106108) 

45. Musacchio A, Salmon ED. 2007 The spindle- 
assembly checkpoint in space and time. Nat. Rev. 
Mol. Cell Biol. 8, 379-393. (doi:10.1038/nrm2163) 

46. Pan J, Chen RH. 2004 Spindle checkpoint regulates 
Cdc20p stability in Saccharomyces cerevisiae. Genes 
Dev. 18, 1439-1451. (doi: 10. 11 01 /gad. 11 84204) 

47. Visconti R, Palazzo L, Grieco D. 2010 Requirement for 
proteolysis in spindle assembly checkpoint silencing. 
Cell Cycle 9, 564-569. (doi : 1 0.41 61 /cc.9.3 . 1 0581 ) 

48. Zeng X et al. 2010 Pharmacologic inhibition of the 
anaphase-promoting complex induces a spindle 
checkpoint-dependent mitotic arrest in the absence 
of spindle damage. Cancer Cell 18, 382-395. 
(doi:1 0.1 01 6/j.ccr.201 0.08.01 0) 

49. Potapova TA, Sivakumar S, Flynn JN, Li R, Gorbsky 
GJ. 2011 Mitotic progression becomes irreversible in 
prometaphase and collapses when Wee1 and Cdc25 
are inhibited. Mol. Biol. Cell 22, 1191-1206. 
(doi:10.1091/mbc.E10-07-0599) 



